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Using a "wall" of quark point sources, we invert the chirally improved Dirac operator to create an "incoherent" 
collection of quark propagators which originate from all spatial points of the source time slice. The lowest-order 
NRQCD approximation is used to create heavy-quark propagators from the same wall source. However, since 
the numerical cost involved in computing such heavy-quark propagators is low, we are able to use a number of 
source gauge paths to establish coherence between the heavy and light quarks at several spatial separations. The 
resulting collection of heavy-light meson correlators is analyzed to extract the corresponding mass spectrum. 



1. Light-quark source and propagators 

For the light quarks, we use the chirally im- 
proved Dirac operator, constructed via an expan- 
sion in gauge paths and the corresponding ap- 
proximate solution of the Ginsparg- Wilson rela- 
tion (for details see Ref. P). 

As sources for the light-quark propagators, we 
use 12 (3 colors x 4 spins) uniform walls, setting 
the color and spin of the quarks to be the same at 
all spatial positions and relying upon the random 
gauges of the ensemble of configurations to cancel 
any non-gauge-invariant (open) contributions to 
hadron correlators constructed from these quark 
propagators. This is similar to the volume source 
used by others [2] for disconnected correlators, 
except for the fact that our source sits on only 
a single time-slice. Connected meson correlators 
from two of these quarks should exhibit a signal 
arising from L 3 closed (or gauge-invariant) cor- 
relators, along with some additional noise from 
the L 3 (L 3 — 1) open correlators. These sources 
have been used before for heavy quarkonium cor- 
relators with some success since the signal-to- 
noise ratio proved to be surprisingly good. Since 
we retain our light-quark propagators and the 
heavy-quark evolution is computationally cheap, 
we are free to use a number of different gauge 
paths to create "walls" of extended heavy-light 



meson sources from the light-quark source before 
we compute the heavy-quark propagators. In this 
way we can construct a full matrix of correlators 
with various spatial extents and a collection of 
starting positions without additional light-quark 
inversions. 

2. Heavy-quark evolution 

The heavy-quark propagators are created using 
lowest-order NRQCD, along with some improve- 
ments for the lattice version 0]. A single time- 
step of heavy-quark evolution is achieved via 
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where we use n — 2 (this is sufficient for b quarks 
on our lattices) and the Hamiltonian is the 0(a 2 )- 
improved, covariant, kinetic-energy operator, cor- 
rected to remove O(a) artifacts in the time evo- 
lution: 
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Since we include only the kinetic-energy term in 
our NRQCD expansion, there is no coupling of 
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the heavy quark's color and spin, the latter being 
provided simply by the projection 

G%. Q (x,0\y,t) = i(l + 7 4 )^G^(f,0|y,t). (3) 

3. Source and sink operators 

The gauge paths and corresponding spin struc- 
tures we use to create our mesons can be found in 
Rcfs. |5I6| . Through sums (s) and differences (p) 
of paths in opposite directions, along with differ- 
ent Dirac matrices, meson operators with differ- 
ent symmetries can be constructed: e.g., 

S : Q(x)TU P (x\x')q(x') =Qj 5 J2 s ^' ^ 

i 

P-'Q^2 HVi 3 ; P+'-Q (HPi - 72P2) q ; (5) 

i 

D±:Q^5 (si - s 2 )q. (6) 

Thus far, we have only included straight paths, 
limiting us to S-, P- and D-wave mesons. Also, 
since we don't include heavy-quark spin interac- 
tions, a number of physical states are not separa- 
ble: the S-wave operator, S, must be viewed as 
creating a spin-average of J p = 0~ and 1~ states; 
the P_, J p = {0+,l+}; the P+, J p = {1+2+}; 
and the D±, J p = {1",2", 3"}. 

4. Heavy-light correlator matrix 

Given all the ingredients, we now present the 
form of our heavy-light correlators: 

Cij(t) = / Yl [T Sa UpW')].G%. q {xAm 

[T%u^(y\y')] j G^(y , ,tr,o)). (?) 

Note the sum over both quark source positions, x 
and x ', due to the use of the walls; only the closed 
terms (where x = x") will contribute to the sig- 
nal. We include both the Qq and Qq correlators 
from opposite time directions. 

After calculating the matrix of correlators, 
Cij (t) , between our linearly independent meson 
operators, tpi = Q(x)TUp(x\x + ri)q(x + ri), we 



solve the corresponding generalized eigenvalue 
problem to extract the physical states [7]: 

]T Cij (t)^ = A fe (t, t ) Cij (to)l# , (8) 

j j 
where, for large enough values of t (3> to), 

\ k (t,t ) = e~ mk{t ~ to) ■ (9) 

So the eigenvalues provide the masses of the 
states and the eigenvectors reveal the meson 
wavefunctions' dependencies upon the Qq separa- 
tion, r{ (= — 7a along the principal axes). This 
technique of using operators with different spatial 
extents has also been recently used to determine 
excited states of baryons |5J . 

5. Preliminary results 

In this section we present results from one set 
(100) of 12 3 x 24 quenched configurations which 
were created using the Liischer-Weisz action at 
13 = 7.9 (r Q /a = 3.384 -> a" 1 » 1330 MeV). 

In Fig. we plot the effective masses calcu- 
lated from the largest eigenvalues of each set of 
operators. The horizontal lines indicate the la 
masses resulting from correlated fits to the data 
over the corresponding time ranges. We repeat 
this for 6 mass combinations (arriQ — 3.0, 3.5 and 
am q = 0.04,0.08,0.10) in order to interpolate to 
the physical values of m b and m s , and to extrap- 
olate to m v u). The physical strange quark mass 
is determined via light-light meson spectroscopy 
{M% k 2M% - Ml -> am s w 0.09); the bottom 
quark mass is determined via finite-momentum 
heavy-light states and the corresponding disper- 
sion relation: 

E(P) « E + ^ , (10) 

where we require that M — (3Mg» + Mg 3 )/4. 
At this lattice spacing we find m b ~ 4240 MeV. 
The interpolated (and extrapolated) results for 
the (fiW) system are displayed in Table 
(Table 13). 

The SP and SD mass splittings are in agree- 
ment with experiment (available only for the for- 
mer) and the findings of previous lattice stud- 
ies |6I1UI11| . All ground-state eigenvector com- 
ponents have the same sign: i.e., there are no 
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Table 1 



Mass splittings for the Bs system 



o 


a(M - M s ) 


M - M s (MeV) 


exp. [5] 


S' 


0.58(11) 


770(150)(?) 




P- 


0.314(28) 


419(37) 


448(16) 


p+ 


0.341(31) 


455(41) 


448(16) 


D± 


0.530(76) 


710(100) 




Table 2 






Mass 


splittings for the system. 




O 


a(M - M s ) 


M - M s (MeV) 


exp. [IP 


5" 


0.46(17) 


610(230)(?) 




P- 


0.306(50) 


408(67) 


384(8) 


P+ 


0.330(58) 


440(77) 


384(8) 


D± 


0.45(11) 


600(150) 





Figure 1. Effective masses versus time from the 
largest eigenvalues of each of the meson operators. 
The results for the P+ and for the first-excited 
state of the S are shifted by 10 in time for clarity. 



0.08, 



3.0 and to /a — 1. 



nodes in the corresponding ground-state heavy- 
light wavefunctions. The first-excited S-wave dis- 
plays a single node at r « 2a w 0.6r (r w 0.5r 
found in 6 ). Our mass results for the first- 
excited S-wave, however, appear to be higher 
than those found previously (« 500 MeV for B s ). 
Our mass plateaus for this state are rather sus- 
pect (typically only 3 t values were used in the 
fits). This issue can only be resolved with bet- 
ter statistics, which can be achieved a number of 
ways: more or larger configurations, more wall 
sources with locally gauge-rotated quark fields, 
more operators (e.g., non-straight paths). We 
would like to consider each of these options in the 
future, along with the inclusion of heavy-quark 
spin interactions to resolve more spin splittings 
in these heavy-light systems. 
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